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THE OVERHAND SHUFFLE MIXES IN @(N 2 LOG AT) STEPS 

By Johan Jonasson 1 
Chalmers University of Technology 

The overhand shuffle is one of the "real" card shuffling methods 
in the sense that some people actually use it to mix a deck of cards. 
A mathematical model was constructed and analyzed by Pemantle 
[J. Theoret. Probab. 2 (1989) 37-49] who showed that the mixing time 
with respect to variation distance is at least of order n 2 and at most of 
order n 2 log n. In this paper we use an extension of a lemma of Wilson 
[Ann. Appl. Probab. 14 (2004) 274-325] to establish a lower bound 
of order n 2 log n, thereby showing that n 2 log n is indeed the correct 
order of the mixing time. It is our hope that the extension of Wilson's 
lemma will prove useful also in other situations; it is demonstrated 
how it may be used to give a simplified proof of the 0(n 3 log n) lower 
bound of Wilson [Electron. Comm. Probab. 8 (2003) 77-85] for the 
Rudvalis shuffle. 

1. Introduction. The distribution of any aperiodic irreducible Markov 
chain approaches the unique stationary distribution as the number of steps 
increases, but how many steps are needed to get sufficiently close? This 
question has attracted much attention over the last two decades or so. One 
important reason for this is that computer development has allowed the de- 
velopment of powerful Markov chain Monte Carlo techniques for simulation. 

Particular interest has been paid to the case when the state space is the 
symmetric group in which case one often refers to the Markov chain as a 
card shuffling chain. A great variety of different card shuffling chains have 
been studied. Among these are many simplified versions of "real" shuffles 
like the top to random shuffle, the transposing neighbors shuffle and the ran- 
dom transpositions shuffle, but also models for shuffles that people actually 
use for shuffling a deck of cards; the riffle shuffle and the overhand shuffle. 
The riffle shuffle (or the Gilbert-Shannon-Reeds shuffle) is a model for the 
most frequently used shuffling method. Aldous [1] showed that the cutoff 
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for a deck of n cards is given by § log 2 n. Later Bayer and Diaconis [2] in a 
celebrated paper extended this to a very sharp result which, among other 
things, established that one needs about seven shuffles to mix an ordinary 
deck with 52 cards. The overhand shuffle, the topic of the present paper, 
is the shuffling technique where you gradually transfer the deck from, say, 
your right hand to you left hand by sliding off small packets from the top of 
the deck. Pemantle [5] introduced and analyzed a probabilistic model for the 
overhand shuffle and established upper and lower bounds of order re 2 log n 
and n 2 , respectively, for the mixing time on a deck of re cards. Here we will 
establish a lower bound of order n 2 logn, thereby proving Pemantle's upper 
bound to be essentially tight. The main tool of our analysis is a technique 
introduced by Wilson [7] where one cleverly uses an eigenvalue close to 1 
and a corresponding eigenvector of the transition matrix of the motion of 
a single card in order to find a useful eigenvector of the whole chain. This 
eigenvector exhibits a certain contraction property and the random variable 
obtained by evaluating it at the state at time t can often be proved to be 
close to its expectation with high probability. In [7] this technique is used 
to prove an order n 3 logn lower bound for the transposing neighbors shuffle 
and in [6] a version for a lifted Markov chain is used to establish a lower 
bound of the same order for the Rudvalis shuffle. Mossel, Peres and Sinclair 
[4] use the technique to prove an re log n lower bound for the so called cyclic 
to random shuffle. 

In Section 2 we review some basic definitions and describe Pemantle's 
model for the overhand shuffle and Wilson's technique. In Section 3 we state 
and prove the promised lower bound for the overhand shuffle. In order to do 
this, we will have to extend Wilson's technique slightly and use a vector that 
is not really but only very close to an eigenvector of the transition matrix. 
We hope that this extension turns out to be useful in other situations as 
well. In Section 4 we demonstrate how it can be used to give a simplified 
proof of Wilson's [6] J7(n 3 logn) lower bound for the Rudvalis shuffle. 

2. Preliminaries. 

2.1. Basic definitions. The standard way of measuring the distance be- 
tween two probability measures \i and v on a finite set S is by the total 
variation norm defined by 

llA* - HI = \ E M S ) " = max(u(A) - v(A)). 

Let {X^}^ , n = 1, 2, 3, . . . , be a sequence of aperiodic irreducible Markov 
chains on state spaces S n and with stationary distributions n n . Assume that 
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\S n \ | oo in a natural way. A sequence {f(n)}^ =1 is said to be an upper bound 
for the mixing time of the sequence of Markov chains if 

]jm||P(l3f (B) e.)-^|=0 

and a lower bound if 

^jnXf (n) E-)-^\\ = l. 

When, for any a > 0, (1 — a)f(n) is a lower bound and (1 + a)f(n) is an 
upper bound, one says that f(n) is a cutoff or a threshold for the mixing 
time. In card shuffling situations, as in the present paper, S n = S n , the set of 
permutations of n cards, and usually 7r n is U = U n , the uniform distribution 
on S n . 

2.2. The overhand shuffle. Pemantle's model for the overhand shuffle is 
parameterized by a probability p G (0, 1). The transition rule is the following: 
Each of the n — 1 slots between adjacent cards is, independently of the 
other slots, declared a outpoint with probability p. This divides in a natural 
way the deck into subsequences, or packets, of cards. Reverse each of these 
packets without changing its position relative to the other packets. (This 
is not what really happens when doing an overhand shuffle; what we get is 
rather the reversed deck compared to what one would have got from the 
overhand shuffle. This of course does not matter for the mixing time.) 

For this Markov chain, Pemantle obtained an 0(n 2 log n) upper bound by 
constructing a coupling for which he was able to demonstrate that one can 
find a constant C < oo such that each card is with high probability matched 
within Cn 2 \ogn steps. Intuitively, one can understand this by observing 
that a particular card makes a random walk which is essentially symmetric 
[at least when the card is not at O(l) distance from the ends of the deck] 
with step size variance of constant order. Therefore, one should be able to 
match the typical card within 0(n 2 ) steps and all cards within 0(n 2 log n) 
steps. Studying the random walks of single cards also quickly yields an Q(n 2 ) 
lower bound since before this time too many of the cards are too close to 
their starting positions. 

2.3. Wilson's technique. The key to the technique is the following lemma 
from [7]. Later on, an extension of the lemma will be needed together with 
corresponding adjustments of its proof. We prove the original lemma here 
in order to make the later proof of the extension easier to digest; doing both 
at once may hide the essential ideas behind technical details. 

Lemma 2.1 (The lower bound lemma). Let {Xt}^ be an aperiodic 
irreducible Markov chain with state space S and stationary distribution ir. 
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Assume that $ : S — > M and 7 G (0, \) are such that 

E[$(X t+1 )|X*] = (l- 7 )$(X t ). 

f7.e., 1 — 7 is an eigenvalue for the transition matrix and & is a correspond- 
ing eigenvector.) Suppose that R > max sg sE [(<!>( A^ + i) — <&(Xt)) 2 \Xt = s], 
put $ = max sg s <]?(s) and let the chain start from a state so such that 

$0o) = 

Fix e > and pui 

r _ log£-(l/2)log(4fl/( 7£ )) 
-log(l-7) 

Then for t <T, 

\\P(X t €-)-ir\\>l-e. 

The lemma may be adopted to situations with complex- valued eigenvalues 
and eigenvectors as well, but we will not need that here; the overhand shuffle 
clearly defines a reversible Markov chain. 

Proof of Lemma 3.2. By induction, E®(X t ) = (1 — 7 )*$(s ) and, con- 
sequently, E^Xoo) = 0. Put A$ = $(X t+ i) - $(X t ). Then 

E[$(X t+1 ) 2 \X t ] = (1 - 2 1 )<$>{X t ) 2 + E[(A$) 2 |X 4 ] < (1 - 2 1 )<$>{X t ) 2 + fl. 

Using induction and that 7 < 1/2, we get 

E[$(^) 2 ]<(1-2 7 )^( S0 ) 2 + ^. 

2 7 

Thus, since 7 < 1/2, 

Var $(X t ) < $( S0 ) 2 ((1 - 2 7 )* - (1 - 7) 2< ) + ^ < ^- 



By Chebyshev's inequality, 



p[ \*(x t )-m(x t )\> J- ) < e - 



7£ 



and by letting t — > 00, it also follows that 



PI |*(*oo)|>l/ — ) < £ 



2 



N ow if f is such that E$(X t ) > 2- v /J?/( 7 e), then we get that P($pf t ) < 
v / P/( 7 £) ) < £/2 and so 

||P(X t G-)-P(^ooG-)ll>l-e, 
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as desired. It remains to determine for which t we have E<3?(A^) = (1 — 7)*$ > 
Some simple algebraic manipulations show that this is the case 
when t <T. □ 

There is more to using the lower bound lemma than just the lemma 
itself: In order to derive good lower bounds, one needs to find an eigenvalue 
sufficiently close to 1 and a useful corresponding eigenvector. It is not at all 
obvious how to do that. In [7], and also in [6] and [4], one may find some 
useful hints on how this may be done. 

3. Lower bound for the overhand shuffle. Let Xt denote the state of the 
deck at time t and let Zt = Z\ denote the position of card i at time t. For 
simplicity, we shall do the the case p = 1/2 and leave the straightforward 
generalization to the reader. To further simplify our calculations, we will to 
begin with assuming a circular deck convention, that is, regarding the top 
card and the bottom card as being next to each other. This means that the 
top packet and the bottom packet may be treated as a single packet (if the 
slot between that top and the bottom card happens to be a cutpoint). Then 
any single card performs a symmetric random walk where the probability 
that the card moves k steps to the right (modulo n) at a single shuffle is 
3(5)^', fcGZ. (A negative k, of course, means that the card moves to the 
left. The given step size distributions, not entirely correct because of the 
remote possibility that, at a given shuffle, there are no cutpoints at all or 
only one cutpoint. However, the probability that this ever happens in the 
time scale we are interested in is so small that we do not need to care about 
this. If we want to make the formula completely correct, we may just redefine 
the shuffle on this event in a suitable way.) We start with the circular deck 
convention since it produces a much cleaner proof than the "real" overhand 
shuffle. At the end of this section the circular convention will be removed 
via a slight extension of the lower bound lemma. 

To find an eigenvector for the {Zt}-chain, we use the cosine function [note 
that cos(2irk/n) is well defined on Z n ]: 



where 7 = (1 + o(l))87r 2 n 2 , which follows from the fact that, for k = o(n), 
one has cos(2vrA;/n) = 1 - 2-K 2 k 2 n~ 2 + 0(k 4 n- 4 ), and that J2T=i k 2 2~ k = 6. 



E cos 



2irZ t+ i 



n 





G 
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The second equality follows from the standard trigonometric formula for the 
cosine of a sum. 

Now put m= \ n/2\ and 

m 9-7T 7 ' 1 

— n 

i=i 

Then <E> is an eigenvector corresponding to the eigenvalue 1 — 7 for the 
transition matrix of the whole deck, and = Cn for C = 0(1). In order to 
use the lower bound lemma, we need to bound K{(Q(Xt+i) — Q(Xt)) 2 \Xt\. 
Write 

27rZ| +1 2ttZ\ 

Yi = cos — - cos -. 

n n 

Then 

m. \ 2 



E[(<t>(X t+1 )-<Z>(X t )) 2 \X t }=E 



E^ 



L \i=l 



X 



Let A denote the event that no two cutpoints of the (t + l)st shuffle are 
more than 10 log n steps apart. It is easily verified that P(A) = 1 — o(n~ 5 ). 

Consider a given pair of cards that at time t are more than 20 log n 
apart. Ideally, if card i and card j were infinitely far apart at time t, then 
Y{ and Yj would have been independent. Put Y( and Y- for such idealized 
independent random variables. By the definition of A, we can couple Yi and 
Yj with Y- and Y- in such a way that Yi = Y( and Yj = Y- on A. (Such 
a coupling can be made by coupling the cutpoints for the two cases in a 
natural way.) Thus, 

ElYiYjlXt] = E\YfY!l A \X t ] + o(n~ 5 ) 



E[Yl\X t ]E[Y^\X t }+o{n 



= E[Yi\X t ]E[Yj\X t }+o{n~ 5 ). 

By mimicking the above calculations, it follows that E[yi|Xj] = 0(n~ 2 ) and 
so 

E[YiY 3 \X t ] = 0{n~ A ). 

Using this, it follows that 

E[{^{X t+l )-^{X t )f\X t ) 

m 

= E E E[K t Y,-|X t ]+0(n- 2 ) 

i=1 j: \Z l t -Z J t \<20logn 
m 

= E E Cov(Yi,Yj\X t ) + O(n~ 2 )<20Vnlogn, 

i=1 j: \Zi-Z{\<201ogn 
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where 



V = maxVar(yi|X f = s) < maxEflf |X t = s] 



8 ^2, sin 2irk/n 9 

for large enough n. Hence, 

E[($(X t+1 ) - <5>{X t )f\X t ] < 2000n~ 1 logn 

and we may apply the lower bound lemma with R = 2000n~ 1 logn. In order 
to have e — > 0, let e = (logn) -1 . The lower bound lemma now yields the 
lower bound 

n 2 / 1 8000?i(logn) 2 



r=(l + o(l))— 5 log(Cn)--log 



8vr 2 V V ' 2 6 8^ 2 
= (l + o(l))^2n 2 logn. 

Finally, for general p, the analogous calculations using that Y^kLi 3^0- ~pV = 
(2 —p)(l —p)/p 3 give the following: 

Theorem 3.1. A lower bound for the overhand shuffle with circular 
deck convention and with parameter p £ (0, 1) is given by 

P 2 (2-P) 2l 



8tt 2 (1-p 2 ) 



n logn. 



Next we remove the circular deck convention. We will again concentrate on 
the case p = 1/2. The vector above is now no longer really an eigenvector 
for {X t }, but only very close to one; in particular, the cards near the ends of 
the deck do not behave well. However, the feeling is that the small error we 
make if we still assume that <I> is an eigenvector should not be large enough 
to upset things too much. This is made precise by the following extension 
of the lower bound lemma: 

Lemma 3.2 (The extended lower bound lemma). Let the setting be as 
in the lower bound lemma with the exception that <I> is not an eigenvector 
for the transition matrix, but there is a positive number p such that 

(1 - -p< E[<t>(X t+1 )\X t ] < (1 - 7)*M + P 

a.s. Then, for t <T, \\P(Xt •) — 7r|| > 1 — where 



log $ - log(2p/ 7 + yJA{R + 6p$)/(7e) ; 
-log(l -7) 
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Proof. The proof is basically a repetition of the proof of the lower 
bound lemma, with the proper adjustments: By induction, we get 

(1 - 7)**(so) - - < E$(X t ) < (1 - 7)*$(*o) + 
7 7 



that is, for some a € [—1,1] 

E$(X t ) = (l-7)**(*o)-- 

7 



and by letting t — > oo, we also get 

-^<E$(X 00 )<^. 

7 7 

We have 

E[<S>(X t ) 2 \X t ] < <5>{X t f + 2$(X t )E[A$|X t ] + R 
<{l-2 1 )<$>{X t ) 2 + 2p\<$>{X t )\ + R 
< (l-2 7 )$(X t ) 2 + 2p$ + R. 

By induction, 

E[$(^) 2 ]<(l-2 7 )^( S o) 2 + ^^. 

Thus, 

Var *(X t ) < (1 - 27) i d>( So ) 2 + - ((1 - 7)^0) - y )' 

^R + 2p$ 2|a|p4>(l -7)* fl + 6p$ 
— 27 7 — 27 

By Chebyshev's inequality, 



p(l* W - E s ( x ( )l>^±Mj<£ 

and by letting t — > 00, using that E$(X 00 ) < p/7, 



Thus, ||P(X t G •) - 71- 1| > 1 - e when t is such that E$(X<) > p/7 + 



2y/(R + 6p$)/('ye). Since E$(X t ) > (1 - 7)*$ -p/7, this holds when t < T. 
□ 
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In our case 7 = 0(n~ 2 ), R = 0(n~ 1 \ogn) (the above arguments apply 
equally well to the present setting without the circular deck convention) and 
<!> = 0(n). Hence, once it has been shown that p = 0(n~ 2+s ) for any 5 > 0, 
then we get the same lower bound that we got with the circular deck con- 
vention. In fact, with some care it can be shown that p = 0(n~ 2 ), but things 
get considerably easier if we settle for proving that p = 0(n _2 (logn) 3 ). 

As above, the probability that in a given shuffle the position of the first 
cutpoint is larger than 10 log n or that the position of the last cutpoint is 
less than n — 10 log n is o(n -5 ). Therefore, when Xt is such that a given card 
is in a position Z t = k with 10 log n < k < n — 10 log n, the movement of that 
card can be coupled with the movement of a card under the circular deck 
convention so that the two movements coincide with probability 1 — o(n -5 ). 
Thus, 



E 



2vrZ 



cos ■ 



t+l 



n 



(1 — 7) cos 



2nZ f 



n 



+ o(n 



Now assume that k < 10 log n (the case k > n - 
same way). Then 

(logn) 2 



10 log n is treated in the 



2-rrk 



cos ■ 



1-0 



Since P(\Z t+ i — Z t \ > 21og 2 n) <n 2 , we get 



E 



2vrZ 



cos ■ 



t+l 



n 



> 



1 



cos ■ 



27r(A; + 21og2n) 



n 



n- 



l-O 



(logn)'' 



(1 — 7) cos 



2irZ t 



11 



+ 



(logn)' 



Since 



it follows on summation that 

E[$(X m )|X t ] = (l- 7 )$(X t 



> cos 

f-f n 

i=i 



+ 20 logn O 
(logn) 3 



(log 



n 



n 



n . _ 
+ -o(n 



= (i- 7 )*M + o 

and so p = 0((logn) 3 ?i -2 ) as desired. 

Again, the case with a general p € (0, 1) is completely analogous and we 
arrive at the following: 
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Theorem 3.3. A lower bound for the overhand shuffle with parameter 
p G (0, 1) is given by 

P 2 (2-P) 2l 

-n log??,. 



8tt 2 (1-p 2 ) 

4. The Rudvalis shuffle. The (inverse) Rudvalis shuffle is described by 
the following transition rule: With probability 1/2, move the bottom card 
to the top of the deck and with probability 1/2, move the card next to the 
bottom card to the top of the deck. The model was proposed by Arunas 
Rudvalis as a shuffle with very slow convergence to uniformity In [3] it was 
shown that the Rudvalis shuffle mixes in 0(n 3 log n) steps. For a long time 
the best known lower bound was Q(n s ), until Wilson [6] recently found an 
0(n 3 logn) bound and thereby established the mixing time as G(n 3 logn). 

Applying the lower bound lemma directly to the Rudvalis shuffle is diffi- 
cult, the reason being that the nontrivial eigenvalues of the chain described 
by the motion of a single card are complex and deviate from 1 by Q(n _1 ). 
Wilson gets around this problem by lifting the Rudvalis chain to a larger 
state space by letting time be a part of this larger state space. He then 
extends the lower bound lemma in a way that lets him make conclusions 
about the mixing time of the Rudvalis chain from the behavior of the lifted 
chain. 

Here we shall reprove Wilson's lower bound by applying our extended 
lower bound lemma above to the n-step transition matrix of the Rudvalis 
shuffle. More precisely, let Xt as before denote the state of the deck at time 
t and Z\ the position of card i at time t. Then with 7 = (1 + o(l))A-K 2 n~ 2 

and ®(X t ) = e!=( 2J cos(27rZ|/n), the pair (l- 7 ,$pf t )) is sufficiently close 
to an eigenvalue/eigenvector pair for the chain {X ns } ( ^ 1 for the extended 
lower bound lemma to be invoked. Due to the fact that we are content with 
finding a vector which is not really but only close enough to an eigenvector, 
we are able to avoid some technical difficulties. For example, we will not 
need to lift the chain and we will not have to deal with any complex algebra 
whatsoever. 

Let pkj denote the probability that a card starting from position k finds 
itself in position j after n steps of the Rudvalis shuffle. If k ^ {n — 1, n}, then 
the first n — k — 1 steps deterministically take the card to position n — 1. 
From there it takes the card a geometric(l/2) number G, say, of steps for 
the card to move to position 1. After that, provided it happens before all 
the n steps have been taken, the remaining k + 1 — G steps deterministically 
shift the card to position k + 2 — G. We get 

Pk,n = (2) 
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and for j = 1, 2, . . . , k + 1, 



A special treatment along the same lines of the case k € {n — l,n} gives 

l i 1 1 \ n 

Pk,n =Pk,l = 4 + (2) 

and for j = 2, 3, . . . ,n — 1, 

Pk }j = (h) n+1 - j - 

The analysis is now very similar to that of the overhand shuffle above. Sup- 
pose X ns is such that Z nt = Z % nt = k for a given card i. When k € {n, n — 1} 
or k < 10 log n, cos(27r£;/n) = 1 — O ((log n) 2 n~ 2 ). Since -P(|Z, 
21og 2 n) < n~ 2 , 



n(s+l) ^ns 



E 



2ttZ, 



cos ■ 



n(s+l) 



7? 



1-0 



(logn) 



(1 — 7) cos 



2ttZ„ 



n 



+ 



(logn) 



11- 



Now assume that k € { [10 log n] , . . . , n — 2} . Then 



E 



2ttZ, 



cos ■ 



n(s+l) 



n 



1 



zZ 92-i 



1 2vr(A;+j) 1 



cos ■ 



+ 



n 



j=-k 

(E t ^-- 

+ o(n~ 5 ). 



2 fc+2 
2vr/c 



cos ■ 



n 



— - sin 



27TJ 



2vr/c 



8111 ■ 



n 



The first of the three terms equals (1 — 7) cos(27r/c/ra). For the second term, 

El . 2irj J^, 1 . 2nj 
-T5 — r sm = > „ , . sm 

j=-k j=-l 



j=_i 



2?r _j_ 

j=-i 



+ 0(n- 3 ) = 0(n" 3 ), 



where the last equality follows from the fact that the last sum is 0, which 
can, for example, be seen by noting that the sum coincides with —2 plus 
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the expectation of a geometric(l/2) random variable. On summing the three 
terms and then over the cards 1,2,..., L n /2j , we get 



EMx n{s+1) )\x ns ] = (i - 7 mx ns ) + o 



(logn) 



In order to apply the extended lower bound lemma, we need to bound 
E[(<S>(X n{s+1) )-<S>(X ns )) 2 \X ns ].Put 

27r(Z l t -t) 



U: = cos ■ 

n 

Then, with m = \n/2\, 

m 

Hx n(s+1) ) - $(x ns ) = Y,(K( S+ i) - Ks) 

i=i 

m n(s+l) 

=E E (ui-uu 

i=l t=ns+l 
n(s+l) m 

= E E(^-^-i)- 

t=ns+l i=l 

Since the terms for different t's are concerned with different shuffles, these 
are independent. Also, since all cards but at most two are at each shuffle 
shifted one step down the deck, at most two terms of the inner sum can 
change, and if so, the change is bounded by 2/n. Taking these facts into 
account, we get 

Var($(X n(s+1) ) - <5>(X ns )\X ns ) < m . — <-. 

Adding this to 



(E[$(X n(fl+1) ) - <S>(X ns )\X t ]) 2 = [ ^(X ns ) + 0[ ) ) = 0(n 



(logn) 



we get 



E[($(JT n(fl+1) ) - <S>(X ns )) 2 \X ns ] = 0{n~ l ). 



We may thus apply the extended lower bound lemma with 7 = (l + o(l))47r 2 ?i -2 , 
R = 0(n _1 ) and p = 0((log?i) 3 n~ 2 ). Doing so, we arrive at the lower bound 
gi^^logn for {X ns }'^ =1 . For the Rudvalis shuffle, this gives the lower bound 



1 

as desired. 



— r^n 3 logn, 



OVERHAND SHUFFLE 



13 



Acknowledgment. I am grateful to the referee for the thorough reading 
of the manuscript. 



REFERENCES 

[1] Aldous, D. (1983). Random walks on finite groups and rapidly mixing Markov 
chains. Seminaire de Probabilites XVII. Lecture Notes in Math. 986 243-297. 
Springer, Berlin. MR770418 

[2] Bayer, D. and Diaconis, P. (1992). Trailing the dovetail shuffle to its lair. Ann. 
Appl. Probab. 2 294-313. MR1161056 

[3] Hildebrand, M. (1990). Rates of convergence of some random processes on finite 
groups. Ph.D. thesis, Harvard Univ. 

[4] Mossel, E., Peres, Y. and Sinclair, A. (2004). Shuffling by semi-random trans- 
positions. FOCS 2001 

[5] Pemantle, R. (1989). Randomization time for the overhand shuffle. J. Theoret. 
Probab. 2 37-49. MR981762 

[6] Wilson, D. B. (2003). Mixing time of the Rudvalis shuffle. Electron. Comm. Probab. 
8 77-85. MR1987096 

[7] Wilson, D. B. (2004). Mixing times of lozenge tiling and card shuffling Markov 
chains. Ann. Appl. Probab. 14 274-325. MR2023023 

Department of Mathematics 
Chalmers University of Technology 
S-412 96 Goteborg 
Sweden 

e-mail : jonasson@math. chalmers . se 



